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DECAY OF HOMOGENEOUS TURBULENCE FROM A GIVEN INITIAL STATE 


by Robert G. Deissler 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 

ABSTRACT 

The homogeneous turbulence problem is formulated by first speci- 
fying the multipoint velocity correlations or their spectral equiva- 
lents at an initial time. Those quantities, together with the corre- 
lation or spectral equations are then used to calculate initial time 
derivatives of correlations or spectra. The derivatives in turn are 
used in time series to calculate the evolution of turbulence quantities 
with time. When the problem is treated in this way the correlation 
equations are closed by the initial specification of the turbulence 
and no closure assumption is necessary. An exponential series which 
is an iterative solution of the Navier-Stokes equations gave much 
better results than a Taylor power series when used with the limited 
available initial data. In general, the agreement between theory and 
experiment was good. 


INTRODUCTION 

A basic difficulty in the usual analyses of homogeneous turbulence 
is the closure problem; that is, the set of correlation or moment equa- 
tions contains more unknowns than equations. The problem occurs, of 
course, because of the nonlinearity of the Navier-Stokes equations 
from which the correlation equations are obtained . 1 
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Although many approximations have been introduced into the corre- 
lation equations (or equivalent spectral equations) in attempts to 
obtain closure, those suggestions have varying degrees of arbitrariness. 
The analyses in Refs. 2 and 3 , although based on definite physical 
ideas, contain dimensionless constants which must be determined by 
experiment. Those in Refs. 4 to 7, although somewhat more deductive, 
have other difficulties. That in Ref. 4, at least for the restricted 

initial condition for which it has been worked out, sometimes gives 

8 

negative energies. The analyses in Refs. 5 and 6 give reasonable 
results for moderately weak turbulence but become unduly complex for 
high Reynolds numbers. That in Ref. 7, although it has yielded some 
realistic deductions, also has computational difficulties because of 
its complexity. 

There is another way of looking at the problem of homogeneous 
turbulence. In order not to lose sight of our goal we will first give 
a statement of that problem. The statement given by Batchelor is 
essentially this; given the statistical state of a homogeneous turbu- 
lent field at an initial instant, to predict the evolution of the tur- 
bulence (in probability) as a function of time. In order to completely 
specify a turbulent field at an initial time, it is necessary to give 
all of the multipoint velocity correlations or their spectral equi- 
valents at that time. 1 It is not hard to show that, given these multi- 
point correlations and the correlation equations, all of the time de- 
rivatives of the turbulent energy tensor and of other pertinent turbu- 
lence quantities can be calculated. These time derivatives can then be 
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used in a series, for instance a Taylor series, to calculate the 
evolution of the turbulent energy tensor (or of the equivalent energy 
spectrum tensor) and of other turbulence quantities. 

It is noted that when the turbulence is treated in this way, there 
is no longer a closure problem. The correlation equations are used 
only to relate the correlations at an initial time to their time de- 
rivatives, and those correlations must be given in order to have a 
complete specification of the turbulence at that time. Of course, in 
practice only a small number of the correlations, and thus of their time 
derivatives will ordinarily be available, but a sufficient number may 
be known to give a reasonably good representation. It might be pointed 
out that even in those analyses which require a closure assumption, the 
turbulence should be specified initially by its correlations or spectra 

since the correlation equations require initial conditions. 
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Kraichnan has very recently studied the convergence properties of 
series such as those considered here. As mentioned in another paper by 
that author' 1 ' 0 , it is not necessary that an expansion be convergent in 
order to be useful, since divergent series can provide excellent asymp- 
totic approximations. 11 

In the present paper we will not concern ourselves primarily with 
convergence questions but will use as a test the agreement with exper- 
iment of the results. Although a Taylor series would no doubt give 
good results if sufficient statistical information were available at 

the initial time, it will be seen that an exponential series which arises 

12 

in a study of the nonlinear decay of a disturbance in a fluid is much 
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more satisfactory when a limited amount of initial information is avail 
able. This is not surprising since the exponential series is an iter- 
ative solution of the Navier-Stok.es equations and thus contains informa 
tion which is not contained in the Taylor series. The resulting formu- 
lation gives results which are in quite good agreement with the avail- 
able experimental data. 


INITIAL TIME DERIVATIVES AND SIMPLE EXPANSIONS 
As mentioned in the introduction, if the multipoint correlations 
are known at an initial instant, as they must be for a complete spec- 
ification of the turbulence at that instant, then the time derivatives 
of the correlations can be calculated from the correlation equations. 
For illustrative purposes we will consider the derivatives of the turbu 
lent energy tensor u^, where u^ and u^ are respectively 
velocity components at the points P and P' separated by the vector 
r, and the overbar indicates an averaged value. Then the first time 

derivative of u.u'. at t = t 1 is given directly by the two-point 
i j x 

5 

correlation equations evaluated at t = t^: 
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r 5u.u'. \ 
i .1 j 

S, ^ J t=t i " ^ 


( u i u j u k)t=t, ( U i u j\)t=t, 
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5FT fiJpVfe. “ 3r7 topt-t. 
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( 1 ) 
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where the pressure-velocity correlations are given by 

1 32 (i^pt =tl 

p &r k S k * d V r i 


( 2 ) 


and a similar equation for 
Eq. (2) is 1 


( V >t=< 


The pertinent solution of 



..A I T ^ \ 

where u'! is the velocity at the point x" = x' - s, and the integrat- 

cj 

ion is over all Is space. This solution is for an infinite fluid, for 

which case the boundary conditions are that pul is bounded for 

r = 0 and zero for r = ^». The quantity p is the density, v is 

the kinematic viscosity and p is the pressure. A repeated subscript 

in a term indicates a summation, with the subscript successively taking 

on the values 1, 2, and 3* The correlation equations are, of course, 

derived from the Navier-Stokes equations. The quantity SuTu lA at 

t=t, can be calculated from Eqs. (l) and (2) if u.u’. and the two- 
r i j 

point triple correlations are known at t=t, . 


The second time derivative of u^u^ is obtained by differentiating 
the two=point correlation equations and evaluating the result at t^. 
This gives 
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The quantity [(d/&t)(u.jU^u£.) in Eq. (3) is obtained from the 

three-point correlation equations'* written for t=t- L and ~r' = 

(The vector "r 1 separates the points P and P".) Thus 



7 



is given by 
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Similar equations are obtained for the other pressure-velocity corre- 
lations. The boundary conditions for Eq. (6) are similar to those for 
Eq. (2). that is, pu*.u£ is bounded for or r^ = 0 and zero for 

in Eq. (3) 

is obtained by letting r • = 0 instead of in Eq. (5). Thus, 

if the turbulence is specified sufficiently well at t=t^ that the 
double, triple and quadruple velocity correlations are known, 

o i 2 

(8 u.u’/ot ) can be calculated. Similarly higher-order derivatives 

1 3 / t=t^ 

are obtained by considering four or more point correlations in the turbu- 
lent field^. With the time derivatives of u^uT” known at t=t^, a 
Taylor series gives u^ui as a function of time as 


Also, an expression for Qa/St)(u u’<ol 
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A similar analysis can be carried out in wave number space. For 
instance, the energy spectrum function E, which shows the contribu- 
tions at various wave numbers to u^u^/2, can be written as 

E = Wfrtj. + (if) t=t x (t " h* + §T (t r V 2 + • • • 

where Be/H is obtained from the Fourier transform of the two-point 
correlation equation (Eq. (9) in ref. 5) &s 

f l(- + [*m^ - 

A 

where dA is an element of surface area of a sphere of radius K> 

K is the wave number vector corresponding to the spacial vector r, 
and and ' are respectively the Fourier transforms of u^u! 

and • Extracting from the integral that portion which can be writ- 

ten in terms of E, and setting the rest of the integral equal to T, gives 
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gf = T - 2 v/c 2 E (10) 

Equation (10) is the well known scalar form of the two-point spectral 
equation. The transfer term T produces energy transfer between wave 
numbers and arises from the triple correlation term in Eq. (l) 

(with i= t j).' 1 ' (Note that the pressure -velocity correlation terms in 
Eq. (l) drop out for i=j.) The second time derivative of E is 




- 2v/C 2 (T) t=t + (2v/C 2 ) 2 (E) += . 


t=t. 


, can be calculated from the two- and three- 

h 

point spectral equations if the two- and three-point spectral quantities 
in those equations are known at t=t^. From Eqs. (20, (23), and (24) 
in reference 5 we obtain 



The quantity (Sy^t)^ 


f ^ijk,^ijki 




dx 1 dA(x) 


—7 n — 

where K is the wave number vector corresponding to r* , 6K = 
dXj dJ <2 d/c 3 , and and are respectively the Fourier 
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transforms of and . If by analogy with the 

procedure used for obtaining Eq. (10 ) , we extract from the integral 
that portion which can, be written in terms of spectral quantities 
already defined (E and T), we have' 


c>T „ 2 m 
+ 2VK T - 


li; 




~ V ( p ijk’ P ijkP 


( 11 ) 


where V is a quantity related to the three-point spectra}, tensors 

^ijk an<i ^i jki ’ More Precisely we can say that V is a functional 

of and since each value of V depends on values of 

^ijk ami ^ijki 811 P°* rrts K ' space. With Eq. (ll), the 

2 A 2 

expression for (5 e/H becomes 




i 


/t*t. 


(V) t=t - 4v* (T) 


t=t. 


+ (2v/C 2 ) 2 (E), . 


t=t. 


( 12 ) 


The Taylor series for E then becomes 


E = + 


(T) t=t - 2 V k (E) t ^ + J (t - t-, ) 
* 1 


:=t j 


ll [j^t^ " 


4v/C 2 (T) t=ti + (2V/C 2 ) 2 (E) t=t 


J 


(t - t .)* 


+ * . 


( 13 ) 
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Equation (13) was used in conjunction with available experimental 
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data at an initial time in an attempt to calculate the variation with 
time of E and thus of uTuT. However, with the available initial 
data ((E) t _ t , (T)^ ^ , and (v)^._ t ), reasonable results were not 
obtained except at small times (Fig. 2). It thus appears that in order 
to obtain good results by using a simple Taylor series, initial 
statistical information of much higher order than that which is avail- 
able would have to be given. Thus, an alternative approach which makes 
more efficient use of the initial statistical information and also 
incorporates additional information from the equations of motion will 
be considered. 

A WORKABLE FORMULATION FOR THE DEVELOPMENT OF 
TURBULENCE FROM A GIVEN INITIAL STATE 
In order to obtain a more efficient means for calculating the 
evolution of turbulence than by a Taylor series in time, we consider 
an iterative solution of the Navi er- Stokes equations similar to that 
in Ref. 12. In addition to the initial statistical information and 
calculated time derivatives we will then have information about the 
form of the decay law from the equations of motion. 

Although attention was confined to determinate initial conditions 
in Ref. 12, for the present purposes we can just as well assume the 
initial velocity fluctuations to be random or turbulent. Thus, we con- 
sider a field of homogeneous turbulence to be made up of a very high 
density of eddies or harmonic disturbances in wave number space. For 
all practical purposes then, since the density of disturbances is very 
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high, the spectrum of the turbulence can be considered continuous. 
The velocity and pressure at any point in the field are given by 


8u. n2 

i . O Ui 

“St “ v 


1 5p a(u iV 

p sr-~s 


*k 


and 

1 5 2 t) _ ^ 

p = - ^ X ; 


(16) 


(17) 


The latter equation is obtained by taking the divergence of Eq. 

( 16 ) and applying the continuity equation. 

From the spectrum of harmonic disturbances we arbitrarily select 
two cosine terms with wave number vectors q and r. Then, the 
velocity associated with those disturbances will be 

u. = a. cos q • • x + b. cos r < x (18) 

where the superscript cc on the velocity indicates that it depends 

on two cosine terms. The results that follow would be the same if two 

sine terms or a sine and a cosine term were considered. If u^ is 

substituted for u^ in the right sides of Eqs, (l6) and (17), the 

time variations of a. and b. plus additional harmonic terms are 

l i 

obtained. If we then substitute that new expression into Eqs. (16) 
and (17), another expression containing still more harmonic terms is 
obtained. In each approximation, the linear terms of the Navier-Stokes 
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equations are considered as unknown and the nonlinear terms as known 
from the preceding approximation. As shown in Ref. 12,. continuation 
of this process leads to 



(19) 


( 20 ) 


( 21 ) 


Comparison of Eqs. (19) to (21) with the first and second approxi- 

c 2 s 2 

mat ions in Ref. 12 shows that t£> ^ = vk and b-*^ = vk . Also, we 
note that since the two harmonic components in Eq. (18) were selected 
arbitrarily, expressions similar to Eqs. (19) to (21) will be obtained 
for any other two components. But the nonlinear interaction of any 
number of harmonic components can be expressed as the sum of the 
interactions of pairs of components (Eqs. (37) and (38), Ref. 12). 

Thus u^, the velocity resulting from all the harmonic components will 
be of the form of equations (19) to (21), and can be written as 
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u. 

1 


D( A iV°° S 

K 


i ■ t 


+ $*«**'*) 


( 22 ) 


where 


A I 

i,7r 


( I „ _ ( ) 4,^-*!) 

/ y i,«»q 
q 

qj*i 


a i,£i e 


(23) 


The summations in Eqs. (22) and (23) will of course contain many orders 

of magnitude more terms than those in Eqs. (19) to (21). Since the 

( ) ( \ 

initial conditions are random, the quantities A, a; A , and 

1 » K ij/C>q 

b4 ■_ are assumed to be random variables. The space -averaged value 
2 

of u^ (no sum on i) is obtained from Eq. (22) by squaring, integrating 
over a cycle, and using the orthogonality property of sines and cosines. 
This gives 


2 

U ±~l_j2 


( 4?) 2 + [<Jj 


(24) 


where 


\i,/c/ V i»K,V 


-2v/C (t-tj 


r.(4, ij 




Lj ' 

q 

q^i 


(25) 
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In obtaining Eq. (25) the various a^ ^ ^ were assumed uncorre- 
lated. According to Eq. (25), <(A? ^) 2 > and <(A? ^.)S> in Eq. (24) 

have the same form, so that we do not need to carry along the super- 
scripts c and s. 

We want to obtain an averaged form of Eq. (25) which is a smoothed 
function of the magnitude of the vector it (but not of its direction) . 
In order to do that, we divide the interval of K = (k.K . )"^ 2 


1 x' 


over 


which disturbances occur into a large number of small increments AC. 


£ 


The terms in U in Eqs. (24) and (25) are divided into groups each of 

K 

which corresponds to a particular Ac. (Note that while the magnitudes 
of the various vectors lying in a particular Ac are approximately 
equal, their directions can of course vary.) The group of terms 
corresponding to each Ac is then subdivided into groups in each of 


which the values of the b. in 

i,/C,q 


£ 

q#L 


do not vary appreciably from 


a value of b (/c). The index s designates a particular increment 

2 i 


in the values of 'the b. 


x,K,q 


Also, for each s. 


average value which we designate by <a. _ > . 

1,/C, » 


i,/C,q 
The summation 


will have an 


in eq. (25), which applies to a particular lt } is then replaced 


£ 

q 

q*L 

by 


XX c i)( a i,?>s (,t) 


-2b s (/c)(t-t 1 ) 


J 
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which applies to a particular AC, and where n 


s,(i) 


is the number 


of terms in which are assigned to the group s for the component 

<1 

q*l 

i. The parenthesis on i indicates that there is no summation on that 

'' 2 

subscript. Then the average value of A . -> in the increment Ac be- 


comes (see Eq. (25)). 




-2v< (t-t 1 ) 




“2b (/c)(t-t 1 ) 
a s 1 (26) 


where n is the number of terms 

rC 


in c 


that lie in AC. 


q 

— 2 17^1 

The expression for u. (Eq. (24)) then becomes 


u i- i) we 


K 


-2vK (t-t^) 


CM) ( a '^V K) 


“2b (/c)(t-t, ) 


(27) 


note that 


To obtain an expression for the energy spectrum function E, we 
1 
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— — 2 2 2 

where u.u. = u n + u n + u_ 

li 1 2 3 

Equations (27) and (28) then give 


( 28 ) 







) 



^VK^t-t^ 

e 


. VA f B s\ (4 ,k s) W 

+ ) , - * ' e 

[_ / \ k/x 


AC 


(29) 


where there is now a summation on i. If Ac is very small, we can 
write, to a very good approximation, 

-2 V K 2 (t-t 1 ) -2b a (#c)(t-t 1 ) 

E(k) = B 2 {k) e + ),B 2 (k:) e (30) 

s 

Equation (30) gives the evolution in time of the energy spectrum func- 
tion from an initial state which is specified by the B's and b’s 
in the equation. 

As shown in the last section, if the turbulence is specified at 
an initial instant, the time derivatives of E can be calculated at 
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that instant by using the Fourier transformed correlation equations. 
Thus, it is desirable to write the B's and b's in Eq. (30) 
in terms of E and its derivatives at the initial time. That can be 
done by evaluating Eq.; (30) and its time derivatives at t = t, 
and solving the resulting system of equations for the B's and b's. 

In what follows, we will first retain only two terms of Eq. 

(30), since we evidently do not have initial experimental data avail- 
able to evaluate the unknown functions in additional terms. (The 
equation resulting from the retention of three terms will be con- 
sidered later.) Equation (30) can then be written conveniently as 


E = (E) 


t=t. 


-2VK (t-t x ) 

C(k) e + (1 


- C) 


-2b(#c)(t-t 1 ) 


(31) 


where 0 < C < 1. 

For C = 1 Eq. (31) reduces to the well-known expression for 
the final period of decay. ^ For the general case (C 4 l) we could 
determine C and b in terms of the first and second derivatives of 
Eq. (31) for t=t, and then evaluate those derivatives by using the 
two-point spectral equations (see Eqs. (10) to (12)). The following 
procedure turns out to be simpler however. By substituting Eq. (31) 
into the spectral Eq. (10) we get for the energy transfer term 


2 -2b(t-t.) 

T = 2(1 - C)(v/C - b)(E), . e 

t-tp 


(T) 


-2b(t-tJ 


l 7 ■ 


t=t. 


e 


(32) 
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Then 


6t 

ST 


2b(T) 


t=t. 



-2b(t-t 1 ) 

e 


( 33 ) 


Comparing the last two members of Eq. (33) and. using Eq. (ll) gives 


(V) t=t 

, .,2 t t l 

b = - 'SBC 


( 34 ) 


From Eqs. (32) and (34) we have 


(T) t=t 

c = 1 ' 


(35) 


Equations (31) and (32) then become 


E = (E) t=t C exp -2v^ 2 (t - t x )] 


+ (1 - C) exp< -2 


V K 


(V) t=t 
2 X t l 




1J 


(t - t. 


( 36 ) 


and 


T = (T) t=t exp 



(V) t=t 

VK 2 - 1 

t=t. 


2(tT 


(■t - V- 


( 37 ) 


From Eq. (11) 
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V = 





(t - t x )| 


( 38 ) 


where C is given by Eq. (35). 

Equations (36) and (37) were obtained by retaining two terms 
on the right side of Eq. (30)* We consider next a higher order 
approximation in which three terms are retained in that equation. 

If Eq. (30), with three terms retained, is substituted into Eq. (10), 
we get for T 


T « 2B ?(K 2 


V 


^(t-t^ 


.+ 2-B^K 2 


-2b (t-t ) 

b 2 ) e 


( 39 ) 


Equation ( 39 ) contains four unknown functions which are to be deter- 
mined by the initial conditions. For that purpose we use Eq. (39) 
and its first three derivatives evaluated at t=t^. Thus we obtain 



( 41 ) 
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B, 


2b 2 T 2 


+ T, 


16b 2 (k 2 - 


b l^ b 2 " b l^ 


(42) 


and 



2b l T 2 * T 3 
l6b|(K 2 - b 2 )(b 1 



( 43 ) 


where T^, T 2 and T^ are the first, second, and third time deriva- 
tives of T at t=t^. The first derivative T-^ can be written in 
terms of the functional (v) t _ t , which gives a representation of 
three point spectral quantities (Eq. (11)). Equations for higher 
order functionals can be obtained by the procedure used for obtaining 
Eq. (ll) for V. Thus by using the four-point spectral equations in 
Ref. 6 (Eqs. (ll) and (12)) we get 

= - 2 v« 2 V + R ( 44 ) 

where R is a functional of three- and four-point spectral quanti- 
ties. Similarly 

|| = - 2 vk 2 r + s (45) 

where S is a functional of three, four, and five-point spectral 
quantities. By using Eqs, (ll), ( 44 ) and ( 45 ), the first, second. 
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and third time derivatives of T at t=t^ in Eqs. (40) to (4-3) can 
be written in terms of higher order spectral quantities as 

h ■ - + 

T 2 - (2vK 2 fm t ^ - + (R) t=ti 

and 

T 3 = -(2v/C 2 ) 3 (T) t=t ^ + 3(2 V/C 2 ) 2 (v) t=t ^ 

- 6 V K 2 (R) t _ t ^ + (S) t-t (48) 

RESULTS AND DISCUSSION 

13 

A comparison between the experimental data of Uberoi and the 
present theory (Eqs. (36) to (38)) is given in Figs. 1 to 4. 14 The 
comparison is made for an initial time corresponding to x/m = 48 in 
the experiment (t* = (v/M 2 )t = 0.00l8l8). (X is the distance down- 
stream from the grid and M is the mesh size of the grid.) For the 
initial specification of the turbulence values of E and T were 
obtained from Figs. 5, 9> and 10 in Ref. 13 . Initial values of V 
were not given directly in Ref. 13 but were estimated from the decay 
data for T and Eq. (ll). Except for experimental error those values 


(^7) 


will be the same as those that might have been measured directly. 

The agreement between the predicted and experimental energy spectra 



24 


for the same initial conditions (Fig. l) appears to be quite good, con- 
sidering the difficulty of the measurements. The calculation of the 
experimental values of E required the differentiation of measured one- 
dimensional spectra and an assumption of isotropy. 

Predicted and experimental values for the decay of mu^ are 
plotted in Fig. 2. . -Thfe 1 agreement between theory and experiment is ' ■ 
excellent for values of t* up to about 0.006. (-Note that spectra 
were measured only for values of t* between 0.00182 and 0.00417). 
Elimination of the moderate deviation for t* > 0.006 might require a 
higher order theory (more terms in Eq. (30) ) , together with additional 
initial statistical information. Alternatively the deviation might be 
due to the amplification at large times of slight inaccuracies in the 
measured initial spectra. The theoretical values for t* less than 
0.00182 were calculated by working backwards from the measured initial 
spectra. Also included in Fig.. 2 is a Taylor series solution which 
uses the same initial information as the exponential series, and the 
curve for the weak turbulence approximation. It might be pointed out 
that the curve for the weak turbulence approximation is not the -5/2 
power decay law usually given for the final period, ^ but is the curve 
obtained by using the measured initial energy spectrum and Eq. (3l) with 
0 = 1 . 


Spectra for the energy transfer term T are plotted in Fig. 3. 

The experimental and theoretical curves are in good agreement except 
near the value of K where (T)^_^. changes sign. The deviation there 
results from a mathematical singularity in Eq. (37) when (T) 


t=t n 


0 . 


However, that deviation does not seem to be serious, because the real 
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physical curve in that region can easily be estimated. This is particu- 
larly true since it is known that the total area enclosed by the T 
spectrum should be zero.''" It appears likely that the difficulty could 
be eliminated if another term were retained in Eq. (30). (More will be 
said about that possibility in the next paragraph. ) The deviation also 


carries over to some extent into the results for E and u^u^. How- 
ever, if one does not use values of K close to the point where (T)^_^. 


changes sign for calculating E and u^u^, "the inaccuracies in those 
quantities will be small. It appears that the overall agreement between 
theory and experiment obtained by using equations (35) to (37) should be 
considered encouraging. 

For the sake of completeness, spectra of the functional V (Eqs. 

(ll) and (38)), the third initial condition specified for the turbulence, 
are plotted in Fig. 4. The agreement between theory and experiment is 
probably within the uncertainty in estimating V from the decay data in 
Ref. 13, except in the vicinity of the point where (T)^^. changes sign. 
Thus the theory predicts the evolution in time of E, T, and V, when 
those quantities are specified at an initial time. 

We have not been able to apply a higher order theory to Uberoi's 
data, that is, to evaluate three instead of two terms in Eq. (30) by 
using the initial data given in his paper. However, we can apply a 
higher order theory to an analysis in Ref. 5, since for that analysis 
we can, in effect, calculate as much initial information as is desired. 
That analysis neglects quadruple correlation terms in the three-point 
correlation equations and should apply, for a particular set of initial 
conditions, at times somewhat before the final period of decay. The 
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initial conditions , as well as values at later times , are given by closed 
form equations in that analysis and thus are better defined than may be 
possible in an experiment. For the present purposes, the analytical 
results from Ref. 5 might in fact be thought of as experimental results 
in which the initial conditions are specified exactly. This is because 
for the model chosen, the analysis is exact, and the initial conditions 
used in both that analysis and the present theory correspond to (or are 
a part of) that model. 

The case considered here corresponds to Fig. 6 of Ref. 5. Values of 
dimensionless E, T, and time derivatives of T for the initial specifica- 
tion of the turbulence (tj = 0.002) are obtained from Eqs. (40) and (39) in 
Ref. 5. We can eliminate the time derivatives of T by introducing V 
(Eq. (ll)) and the higher order functionals R and S (Eqs. (44) and (45)). 
In the present case, those quantities will all be representations of corre- 
lations of order no higher than the third, since terms involving correla- 
tions of higher order than the third are assumed negligible in the analysis 
of Ref. 5. 

Figure 5 gives a comparison between results for T calculated from 
the present analysis and those from Ref. 5. The quantity Jq is a con- 
stant related to conditions at tg = -0.00633 in the equations of Ref. 5. 
The starred quantities in Figs. 5 to 7 are the same as those in Figs. 1 to 4 
if we let Jq = . As expected, when T is calculated from Eq. (37), 

the agreement with Ref. 5 is good except in the region where (T)^._^. 
changes sign. However, when a higher order theory is used by retaining 
three terms in the expression for E (two terms in expression for T) (Eq. 
(39)) the agreement is excellent at all values of K. It might be expected 
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that a similar improvement would be obtained in Figs. 3 and 4 if a higher 
order theory could be used for comparison with the experimental data of 
Uberoi. 

Because of the good agreement obtained for T in Fig. 5, one would 
expect the calculated energy spectra E to also be in good agreement with 
those from Ref. 5. Figure 6 shows that that is indeed the case. In order 
to show the effects of energy transfer between wave numbers, curves for the 
final period of decay (first term of Eq. (40} of Ref. 5) are also included 
in Fig. 6. 

Figures 7 to 9 show plots for the decay of the higher-order spectral 
quantities V, R, and S. The agreement between the present higher order 
theory and the results of Ref. 5 is very good. Thus by specifying the 
initial conditions for E, T, V, R, and S, we can predict the evolution 
in time of those quantities by using the present higher order theory. That 
is, the required number of initial conditions is no greater than the number 
of quantities whose decay we can predict. 

CONCLUDING REMARKS 

If a homogeneous turbulent field is specified at an initial instant by 
its multipoint-velocity correlations (or their spectral equivalents), the 
initial time derivatives of those quantities can be calculated from the 
correlation or spectral equations. The development of the turbulence in 
time can then be obtained by using those derivatives in a series such as 
a Taylor power series. When the problem is formulated in this way, an 
assumption for closing the system of correlation equations is not required, 
since those equations are closed by the initially specified correlations or 
spectral quantities. A Taylor series expansion, however, did not give 
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realistic results (except for small times) when the limited initial ex- 
perimental data were used. An exponential series (Eq. (30)), which is an 
iterative solution of the Navier-Stokes equations worked much better with 
the limited initial information. 

In general, when the energy and transfer spectra and a quantity 
related to three-point spectra were specified at an initial time, the 
predicted changes with time of those spectra, as well as the:, 
turbulent energy, were in good agreement with experiment. Since the pre- 
diction of the changes of those spectra with time is evidently an essen- 
tial part of the homogeneous turbulence problem, the results are encour- 
aging . 

A higher-order theory was given in which the above quantities, as 
well as two additional higher order spectral quantities, were specified 
initially. Very good agreement was obtained between the predicted decay 
of all of those quantities and the results for a previous analytical model 
(Ref. 5) . For the present purposes the results for the previous model 
might be thought of as experimental results in which the initial conditions 
are specified exactly. Thus when the results from the present, theory are 
compared with either experimental results or the results of an "analytical 
experiment," the agreement is good, and the number of specified initial 
conditions need be no greater than the number of quantities whose decay 


we can predict. 
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Figure 1. - Comparison of theory and experiment for de 
cay of turbulent energy spectra. 



Figure 2. - Comparison of theory and experiment for 
decay of average component of velocity variance. 
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Figure 3. - Comparison of theory and experiment for 
decay of energy transfer spectra. 
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Figure 4. - Comparison of theory and experiment for decay 
of higher order spectral quantity V (Eq. (ID). 
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Figure 5. - Comparison of present theory with that of 
reference 5 for decay of energy transfer spectrum. 



PRESENT THEORY 

REFERENCE 5 (QUADRUPOLE 

CORRELATION TERMS 
NEGLECTED) 

FINAL PERIOD SOLUTION 



Figure 6. - Comparison of present 
theory with that of reference 5 
for decay of turbulent energy 
spectrum. 



Figure 7. - Comparison of present (higher 
order) theory with that of reference 5 for 
decay of the higher order spectral quantity 
V (Eq. (ID). 
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Figure 8. - Comparison of present (higher order) 
theory with that of reference 5 for decay of the 
higher order spectral quantity R (Eq. (44)). 
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Figure 9. - Comparison of present (higher order) 
theory with that of reference 5 for decay of the 
higher order spectral quantity S (Eq. (45)). 
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